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TITLE OF THE INVENTION 

^ELECTROMAGNETIC SURVEYING FOR HYDROCARBON RESERVOIRS 

BACKGROUND OF THE INVENTION 

The invention relates to seafloor electromagnetic (EM) surveying for oil and 
5 other hydrocarbon reserves. 

Geophysical methods for mapping subterranean resistivity variations by 
various forms of EM surveying have been in use for many years [1, 2, 3, 10], In these 
methods, electric field detectors are placed on the seafloor at carefully chosen 
positions at ranges up to about 10 km from an electromagnetic source. Detector 
10 signals measured at the detectors are sensitive to variations in subterranean strata 
configuration resistivity beneath the area being surveyed. However, EM surveying 
was not widely thought of as a technique that could be applied to finding hydrocarbon 
reservoirs. 



15 reservoirs. An early proposal by Statoil was to use the vertical current flow 
components to detect hydrocarbon layers [4, 5], since it is these components that are 
sensitive to the presence of a thin resistive layer. This was based on the understanding 
that a subterranean strata configuration that includes a resistive hydrocarbon layer 
embedded within less resistive sediments will give rise to a measurable enhancement 

20 of the electric field amplitude compared to a subterranean strata configuration 
comprising only water-bearing sediments. The Statoil proposal was to collect data 
from detector locations which are in-line with (i.e. end-on to) the axis of a horizontal 
electric dipole (HED) antenna so that the galvanic mode, that should be most sensitive 
to the presence of a buried high resistivity layer, dominates. However, it was 

25 established that the Statoil method could not provide reliable results, since the in-line 
data collected is incapable of distinguishing between a thin buried hydrocarbon layer 
of high resistivity situated in less resistive strata, on the one hand, and a non- 
hydrocarbon bearing rock formation in which the strata exhibits increasing resistivity 



More recently, it was proposed to use EM surveying to find hydrocarbon 
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with depth, on the other hand, the latter being a common feature of many large scale 
sedimentary structures. 

It was then proposed to use the EM surveying method according to Sinha [12] 
for finding hydrocarbon reservoirs [9, 13] and it was then confirmed that this method 
5 works well in practice for finding hydrocarbon reservoirs [6, 7]. The essence of the 
Sinha method is to normalise the in-line data with equivalent data for the same 
source-detector pair locations collected in an orthogonal geometry where the inductive 
mode dominates the response, referred to as a broadside geometry. In the broadside 
geometry the axis of the HED dipole antenna of the source is perpendicular to a line 

10 between the detector and source. EM surveying of a hydrocarbon reservoir applying 
the Sinha method is now described in more detail. 

Survey data is collected by using a surface vessel to tow a submersible vehicle 
carrying a HED antenna over a survey area. The HED antenna broadcasts a source 
electromagnetic signal into the seawater. Detectors are located on the seafloor over the 

15 survey area and measure a signal in response to EM fields induced by the HED 
antenna. The amplitude of the detector signals is sensitive to resistivity variations in 
the underlying strata configuration and this is used to determine the nature of the 
subsea structure. In order to successfully map subterranean resistivity variations, the 
orientation of the current flows induced by the source electromagnetic signal must be 

20 considered [6]. The response of seawater and subterranean strata to the source 
electromagnetic signal is different for horizontally and vertically flowing induced 
current components. For horizontally flowing current components, the coupling 
between the layers comprising the subterranean strata is largely inductive. This means 
the presence of a thin resistive layers (which is indicative of a hydrocarbon layer) does 

25 not significantly affect the detector signal at the seafloor since the large scale current 
flow pattern is not affected by the thin resistive layer. On the other hand, for vertical 
current flow components, the coupling between layers is largely galvanic (i.e. due to 
the direct transfer of charge). In these cases even a thin resistive layer strongly affects 
the detector signals at the seafloor since the large scale current flow pattern is 

30 interrupted by the resistive layer. 
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While it is the vertical components of induced current flow which are most 
sensitive to the presence of a thin resistive layer, sole reliance on these components 
for detecting a hydrocarbon layer is not possible without ambiguity. The effects on the 
amplitude signals at the detectors arising from the presence a thin resistive layer can 
5 be indistinguishable from the effects which arise from other realistic large scale 
subterranean strata configurations. In order to resolve these ambiguities, it is necessary 
to determine the response of the subterranean strata to both horizontal (i.e. inductively 
coupled) and vertical (i.e. vertically coupled) induced current flows [6], 

An electromagnetic source such as a HED antenna generates both inductive 

10 and galvanic current flow modes, with the relative strength of each mode depending 
on source-detector geometry. At detector locations which are broadside to the HED 
antenna dipole axis, the inductive mode dominates the response. At detector locations 
which are in-line with the HED antenna dipole axis, the galvanic mode is stronger [6, 
8, 9, 10]. Accordingly, the response of the subterranean strata to vertical induced 

15 current flows along a line between a source location and a detector location is 
determined by arranging the HED antenna to present an end-on orientation to a 
detector, and the response of the subterranean strata to horizontal induced current 
flows along the line between the source location and the detector location is 
determined by arranging the HED antenna to present a broadside orientation to the 

20 detector. Data from both geometric configurations is required. 

It is therefore important when designing a practical EM survey for detecting 
buried hydrocarbon layers using known techniques to distinguish between source and 
detector configurations in which the coupling between layers is largely inductive due 
to horizontal currents (in which case the survey has little sensitivity to the presence of 

25 a thin resistive layer) and those in which the coupling between layers is largely 
galvanic due to vertical currents (in which case blocking of the passage of this current 
flow by a reservoir leads to a survey which is strongly sensitive to the presence of a 
thin resistive layer). 

Figure 1 shows in plan view an example survey geometry according to the 
30 Sinha method. There are sixteen detectors 25, and these are laid out in a square grid 
on a section of seafloor 6 above a subterranean hydrocarbon reservoir 56. The 
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hydrocarbon reservoir 56 has a boundary indicated by a heavy line 58. The orientation 
of the hydrocarbon reservoir is indicated by the cardinal compass points (marked N, 
E, S and W for North, East, South and West respectively) indicated in the upper right 
of the figure. To perform a survey, a source such as a HED antenna, starts from 
5 location 'A* and is towed along a path indicated by the broken line 60 through 
location 'B* until it reaches location *C\ which marks the end of the survey path. As 
is evident, the tow path first covers four parallel paths aligned with the North-South 
direction to drive over the four "columns" of the detectors. This portion of the survey 
path moves from location 4 A* to location *B\ Starting from location *B\ the survey 

10 path then covers four paths aligned with the East- West direction which drive over the 
four "rows" of detectors. Each detector is thus passed over in two orthogonal 
directions. The survey is completed when the source reaches the location marked 'C\ 

During the towing process, each of the detectors 25 presents several different 
orientation geometries with respect to the source. For example, when the source is 

15 directly above the detector position Dl and on the North-South aligned section of the 
tow path, the detectors at positions D5, D6 and D7 are at different ranges in an end-on 
position, the detectors at positions D2, D3 and D4 are at different ranges in a 
broadside position and the detector at positions D8 and D9 are midway between. 
However, when the source later passes over the detector position Dl when on the 

20 East-West aligned section of the tow path, the detectors at positions D5, D6 and D7 
are now in a broadside position, and the detectors at position D2, D3 and D4 are in an 
end-on position. Thus, in the course of a survey, and in conjunction with the 
positional information of the source, data from the detectors can be used to provide 
details of the source electromagnetic signal transmission through the subterranean 

25 strata for a range of distances and orientations between source and detector. Each 
orientation provides varying galvanic and inductive contributions to the signal 
propagation. In this way the continuous towing of the source can provide a survey 
which samples over the extent of the subterranean reservoir. 

The Sinha method has been demonstrated to provide good results in practice. 

30 However, it has some limitations. 
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Firstly, since the two modes cannot be easily separated there will generally be 
a level of cross-talk between them at a detector and this can lead to ambiguities in the 
results. 

Secondly, in order to obtain survey data from both in-line and broadside 
5 geometries, the HED antenna needs to present two different orientations at each 
source location. This requires the surface vessel to make multiple passes over 
broadcast locations and can lead to long and complex tow path patterns. 

Thirdly, the survey can only provide the best data possible at discrete source 
locations. This is because of the geometric requirements of a HHP antenna survey 

10 which dictate that, at any point during the survey, data can only be optimally collected 
from those detectors to which the HED antenna is arranged either in-line or broadside. 
At other orientations, separation of the inductively and galvanically coupled signals 
becomes more difficult, and resulting data are less reliable. For instance, referring to 
the figure, when the HED antenna is at a point on the tow path directly above the 

15 detector marked Dl and on the North-South aligned section of the tow path, in-line 
data can only be collected from the detectors marked D5, D6 and D7, whilst broadside 
data can only be collected from the detectors marked D2, D3 and D4. The other 
detectors (for example those marked D8, D9 and D10) provide only marginally useful 
information at this point of the survey because of the complex mixing of the 

20 galvanically and inductively coupled modes. Furthermore, if, for example, the HED 
antenna is at the location identified by reference numeral 57 in the figure, which is on 
a North-South aligned section of the tow path, in-line data can be collected from the 
detectors marked D3, D8, D9 and D10, but broadside data cannot be collected from 
any of the detectors. Since both broadside and in-line data are required for optimal 

25 analysis, the best data possible with the square detector array shown in the figure can 
only be collected from points along the tow path where the source is directly above 
one of the detector locations. 

In summary, with the Sinha method, the time during which good quality data 
can be collected represents only a small fraction of the overall time taken to perform a 

30 survey. Furthermore, in addition to the survey being time-inefficient, it is necessary to 
accurately follow a complex tow path which has to complement the detector layout, 
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and the detectors themselves must also be carefully accurately arranged. The 
difficulties in controlling both the position and the orientation of a towed source 
antenna, coupled with this need to accurately follow a particular tow path relative to 
the detector grid, is one of the major sources of error in surveys of these kind. The 
5 disadvantages associated with the survey constraints imposed by the Sinha method are 
the price to pay for resolving the ambiguities inherent in the Statoil method. 
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SUMMARY OF THE INVENTION 

According to the invention there is provided an electromagnetic survey 
method for surveying an area that is thought or is known to contain a subterranean 
5 hydrocarbon reservoir, comprising: transmitting a source electromagnetic signal from 
a source location; detecting a detector signal at a detector location in response thereto; 
and obtaining survey data indicative of phase difference between first and second 
components of the detector signal resolved along first and second directions 
respectively. 

10 By comparing phase measurements of different components of the detector 

signal, a phase separation anomaly can be detected which is sensitive to the presence 
of a hydrocarbon layer or reservoir within a subterranean strata configuration. The 
presence or not a phase separation anomaly, and hence the presence or not of a 
hydrocarbon layer, can be determined with a single source orientation. There is no 

15 need, as there is with known methods based on amplitude, for data to be collected 
with different source orientations. Accordingly, surveys can be performed more 
quickly and without needing to accurately control the source orientation. Furthermore, 
because of this insensitivity of a phase measurement to the relative source orientation, 
reliable data collection is not limited to specific source location and detector location 

20 geometries, as is the case when collecting in-line/broadside amplitude data, and a 
much less complex and geometrically restrained towpath can be employed to survey 
an extended area. 

The first and second components can be any two of radial, vertical and 
azimuthal. The clearest phase anomaly appears to occur from the pairing of radial and 
25 azimuthal components. It is also possible to use all three components together, i.e. to 
have first, second and third components. 

The first, second and, if used, third directions are preferably orthogonal, since 
by observing geometrically independent components of the detector signals, there is 
minimal cross-talk between the first and second data sets, and the sensitivity to the 
30 presence of a hydrocarbon reservoir is accordingly increased. 



WO 2004/049008 



PCT/GB2003/005094 



The source electromagnetic signal can be broadcast from an antenna mounted 
on a submersible vehicle, or from a static location, such as within a borehole, or from 
an oil or gas platform. 

The source electromagnetic signal can be emitted at different frequencies to 
5 obtain survey data at a plurality of different frequencies. Moreover, the source 
electromagnetic signal can be emitted at a variety of frequencies, preferably between 
0.01 Hz and 10 Hz. The method can be advantageously repeated over the same survey 
area using different frequencies of source electromagnetic signal. Lower frequencies 
are generally preferred. By probing the subterranean strata at a number of different 
10 frequencies of source electromagnetic signal, it is possible to obtain improved vertical 
resolution of structures within the subterranean strata configuration. 

The source signal can be from a horizontal electric dipole. Such a signal can be 
provided using existing equipment, and also allows relatively simple inversion 
modelling. 

15 The invention also provides a method of analysing results from an 

electromagnetic survey of an area that is thought or known to contain a subterranean 
hydrocarbon reservoir, comprising: providing survey data indicative of phase 
difference between first and second components of a detector signal resolved along 
first and second directions respectively; extracting the phase differences from the 

20 survey data; and determining a metric from the phase differences that is predictive of 
the presence or absence of hydrocarbon. 

The phase differences can be extracted by rotationally transforming the survey 
data from an instrument frame to a source frame. 

The invention also provides a computer program product bearing machine 

25 readable instructions for implementing the analysis method. 

The invention further provides a method of planning an electromagnetic 
survey of an area that is thought or known to contain a subterranean hydrocarbon 
reservoir, comprising: creating a model of the area to be surveyed including a seafloor, 
a rock formation containing a postulated hydrocarbon reservoir beneath the seafloor, 

30 and a body of water above the seafloor; setting values for depth below the seafloor of 
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the postulated hydrocarbon reservoir and resistivity structure of the rock formation; 
and performing a simulation of an electromagnetic survey in the model to obtain from 
the model phase differences between first and second components of a detector signal 
resolved along first and second directions respectively. 
5 Repeated simulations for a number of distances between a source and a 

detector and frequencies can be performed in order to allow optimum surveying 
conditions in terms of source-to-detector distance and frequency of EM signal for 
probing the hydrocarbon reservoir to be selected when performing an electromagnetic 
survey. The effects of differing detectors array configurations and source tow paths 
10 can also be modelled. 

The invention also provides a computer program product bearing machine 
readable instructions for implementing the planning method. 



15 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of the invention and to show how the same may be 
carried into effect reference is now made by way of example to the accompanying 
5 drawings in which: 

Figure 1 is a schematic plan view showing an example survey geometry 
following prior art principles in which sixteen detectors are laid out on a section of 
seafloor above a subterranean reservoir; 

Figure 2A shows in schematic vertical section a surface vessel undertaking an 
10 EM survey; 

Figure 2B is a plan view detailing a polar coordinate system; 

Figure 3 shows in schematic vertical section a model uniform background 
subterranean strata configuration; 

Figure 4 shows in schematic vertical section a model hydrocarbon-layer 
15 subterranean strata configuration; 

Figure 5A shows a graph plotting calculations of phases of different 
components of detector signals seen during a model electromagnetic survey of the 
subterranean strata configurations shown in Figures 3 and 4; 

Figure 5B shows a graph plotting differences in the phases shown in Figure 

20 5A; 

Figure 6 shows in schematic vertical section a model of a non-hydrocarbon 
bearing subterranean strata configuration; 

Figure 7 shows a graph plotting calculations of phases of different components 
of detector signals seen during a model electromagnetic survey of the subterranean 
25 strata configurations shown in Figures 3, 4 and 7; 

Figure 8 shows a graph plotting differences in the phases shown in Figure 7; 

Figure 9A is a schematic plan view showing an example survey geometry 
according to an embodiment of the present invention in which sixteen detectors are 
laid out on a section of seafloor above a subterranean reservoir; 
30 Figure 9B compares the signal coverage of the prior art method and the 

method of the invention; 
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Figure 10 shows a graph plotting calculations of phases of different 
components of detector signals seen during a model electromagnetic survey of the 
subterranean strata configurations shown in Figures 4 at two different source 
electromagnetic signal frequencies; 
5 Figures 11A and 11B show graphs plotting differences in the phases shown in 

Figure 10; 

Figure 12A shows a graph plotting calculations of differences in phases 
between different components of detector signals seen during a model electromagnetic 
survey of several hydrocarbon-layer subterranean strata configurations with an 
10 electromagnetic source signal frequency of 0.5 Hz; 

Figure 12B shows a graph plotting calculations of differences in phases 
between different components of detector signals seen during a model electromagnetic 
survey of several hydrocarbon-layer subterranean strata configurations with an 
electromagnetic source signal frequency of 0.25 Hz; 
15 Figure 12C shows a graph plotting calculations of differences in phases 

between different components of detector signals seen during a model electromagnetic 
survey of several hydrocarbon-layer subterranean strata configurations with an 
electromagnetic source signal frequency of 1.0 Hz; 

Figure 12D shows a graph plotting calculations of differences in phases 
20 between different components of detector signals seen during a model electromagnetic 
survey of several hydrocarbon-layer subterranean strata configurations with an 
electromagnetic source signal frequency of 2.0 Hz; 

Figure 13 A shows a graph plotting calculations of phases of different 
components of detector signals seen during a model electromagnetic survey of several 
25 uniform background subterranean strata configurations; 

Figure 13B shows a graph plotting differences in the phases shown in Figure 

13A; 

Figure 14 shows a graph plotting calculations of maximum observed 
differences in phases between different components of detector signals seen during a 
30 model electromagnetic survey of a hydrocarbon-layer subterranean strata 
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configuration, as a function of hydrocarbon-layer resistivity, and for several 
electromagnetic source frequencies; 

Figure 15 A shows a graph plotting calculations of phases of different 
components of detector signals seen during a model electromagnetic survey of the 
5 subterranean strata configurations shown in Figures 3 and 4; 

Figure 15B shows a graph plotting differences in the phases shown in Figure 

15 A; 

Figure 16A shows a graph plotting calculations of phases of different 
components of detector signals seen during a model electromagnetic survey of the 
10 subterranean strata configurations shown in Figures 3 and 4; and 

Figure 16B shows a graph plotting differences in the phases shown in Figure 

16A; 



15 
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DET AILED DESCRIPTION 

A method of electromagnetic surveying for oil and other hydrocarbon reserves 
is described which does not require separate data acquisition of the response of a 

5 subterranean strata configuration to inductively and galvanically coupled modes. The 
new method can be performed using pre-existing survey equipment. 

Figure 2A schematically shows a surface vessel 14 undertaking EM surveying 
of a subterranean strata configuration in a way that is suitable for collecting survey 
data for carrying out the invention. The subterranean strata configuration includes an 

10 overburden layer 8, an underburden layer 9 and a hydrocarbon layer (or reservoir) 12. 
The surface vessel 14 floats on the surface 2 of the seawater 4. A deep-towed 
submersible vehicle 19 carrying a HED antenna 21 is attached to the surface vessel 14 
by an umbilical cable 16 providing an electrical, optical and mechanical connection 
between the deep-towed submersible vehicle 19 and the surface vessel 14. The HED 

15 antenna broadcasts a source electromagnetic signal into the seawater 4, 

One or more remote detectors 25 are located on the seafloor 6. Each detector 
25 includes an instrument packages 26, a detector antenna 24, a floatation device 28 
and a ballast weight (not shown). The detector antenna 24 measures a detector signal 
in response to EM fields induced by the HED antenna in the vicinity of the detector 

20 25, the amplitude of the detector signals is sensitive to resistivity variations in the 
underlying strata configuration. The instrument package 26 records the detector 
signals for later analysis. The detector antenna 24 in this example comprises two 
orthogonal dipole antennae arranged to detect first and second components of the 
electric field in a horizontal plane, i.e. one which is parallel to the seafloor 6. 

25 The detectors record two (or three) orthogonal components of the seafloor 

electric field as raw data. The raw survey data are then analysed, after recovery of the 
detectors and transfer of the raw data into a suitable computer. Initially a spectral 
analysis is performed to remove the component of the signal which corresponds to 
source transmission, as is conventional. The survey data are then combined with 

30 source and receiver navigation data, again as is conventional. Then the survey data are 
processed to rotate the electric fields from an 'instrument' frame (i.e. components 



WO 2004/049008 



PCT/GB2003/005094 



-14- 

parallel to the receiver dipoles of the detector) to the 'source' frame (i.e. radial and 
azimuthal components referenced to the source-receiver geometry). This is a new 
processing step specific to the present invention. 

Figure 2B is a schematic plan view detailing a polar-coordinate system which 
5 is used to describe the principles of the new method. The origin of the coordinate 
system is positioned at the centre of the HED antenna shown in Figure 2A, and zero- 
azimuth is aligned parallel to the dipole axis of the HED antenna, as indicated in 
Figure 2B (the HED antenna in this Figure not drawn to scale). In Figure 2B, a single 
detector 25 is shown positioned at a range of R km from the origin, and at an azimuth 

10 of 0°. The orthogonal dipole antennae comprising the detector antenna 24 are 
arbitrarily oriented in the horizontal plane as indicated in the figure. 

Since the phase difference from the source of electromagnetic radiation from a 
horizontal dipole source is an azimuthally symmetric subterranean strata configuration 
is largely independent of azimuth O, phase measurements recorded at the detector 

15 antenna 24 are largely independent of the azimuthal position of the detector 25 shown 
in Figure 2B. This allows phase data to be collected equally over a wider range of 
source-detector orientations than is possible with amplitude data, and any inaccuracies 
in the measurement of the azimuthal position of the detector in the coordinate system 
shown in Figure 2B have a lesser effect. 

20 In the following examples, the two components of detected electric field, also 

known as detector signal, for which the phase is measured are a radial component and 
an azimuthal component. The radial component is that component of the electric field 
resolved along a direction parallel to a line connecting the source location and the 
detector location, and marked E p in Figure 2B. The azimuthal component is that 

25 component of the electric field resolved along a direction perpendicular to a line 
connecting the source location and the detector location and in a horizontal plane, and 
marked Eo in Figure 2B. The components of the detected electric field along these 
directions is determined from the angular orientation of the orthogonal dipole 
antennae comprising the detector antenna 24 relative to the line joining the source 

30 location and the detector location. This can be easily determined using standard 
instrumentation, such as, for example, active or passive sonar to determine the relative 
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positions of the source location and the detector location, and a magnetic compass to 
determine the detector antenna orientation. 

In order to show how the respective phases of two spatial components (e.g. 
radial p and azimuthal © components) of the electric field can be used to detect the 
5 presence of a subterranean hydrocarbon reservoir, numerical forward modelling of the 
kind described by Chave and Cox [11] is applied to different model subterranean 
strata configurations. 

Figure 3 shows in schematic vertical section a model background subterranean 
strata configuration. The configuration comprises a section of seafloor 106 beneath a 

10 10 km depth of seawater 104. The seawater has a resistivity of 0.31 Sim. Beneath the 
seafloor 106 is a uniform half-space sedimentary structure with a resistivity of 1 Sim, 
the low resistivity being primarily due to aqueous saturation of pore spaces. This 
background subterranean strata configuration extends uniformly downwards for an 
infinite extent. Also indicated in Figure 3 are a HED antenna 21, and a detector 25, 

15 such as those shown in Figure 2 A. The distance between the HED antenna and the 
detector (i.e. the range) is R km. The azimuthal position of the detector relative to the 
orientation of the HED antenna is arbitrary due to the insensitivity of the phase 
component of the detected electric field signals to azimuth. 

Figure 4 shows in schematic vertical section a model hydrocarbon-layer 

20 subterranean strata configuration. A section of seafloor 106 lies beneath a 10 km 
depth of seawater 104 which has a resistivity of 0.31 Sim. The strata configuration 
beneath the seafloor 106 comprises a 1 km thick overburden layer 108, representing 
sediments, arranged above a hydrocarbon layer 112. The overburden layer 108 has a 
resistivity of 1 Sim, again, primarily due to aqueous saturation of pore spaces. The 

25 hydrocarbon layer 112 is 0.1 km thick, and has a resistivity of 100 Sim. The relatively 
high resistivity of the hydrocarbon layer is due to the presence of non-conducting 
hydrocarbon within pore spaces. Below the hydrocarbon layer 112 is a sedimentary 
underburden layer 109, which, as for the overburden layer, has a resistivity of 1 Sim. 
The underburden layer extends downwardly for an effectively infinite extent. 

30 Accordingly, except for the presence or absence of the hydrocarbon layer 112, the 
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background subterranean strata configuration of Figure 3 and the hydrocarbon-layer 
subterranean strata configuration of Figure 4 are identical. A HEP antenna 21 and a 
detector 25 are again shown as in Figure 3. 

Figure 5A shows a graph plotting the modelled phase 0 of the radial and 
5 azimuthal components of the detected electric field for both the background 
subterranean strata configuration and the hydrocarbon-layer subterranean strata 
configuration models shown in Figures 3 and 4 respectively as a function of range R. 
The phase is measured relative to a source electromagnetic signal transmitted by the 
HED antenna 21. In this example, the source electromagnetic signal is at a frequency 
10 of 0.5 Hz, The radial and azimuthal components of the detected electric field for the 
background subterranean strata configuration are marked 0 P B and 0o B respectively and 
the corresponding components of the detected electric field for the hydrocarbon-layer 
subterranean strata configuration are marked 0 P R and 0<& R respectively. The results 
show that 9 P B , 0<d B , 0 p r and 0<r> R all advance steadily in phase with increasing range R. 

15 However, it is also clear that the rate of phase advance is less for the hydrocarbon- 
layer subterranean strata configuration than for the background subterranean strata 
configuration. In the case of the background subterranean strata configuration, the 
phase of both the radial and azimuthal components advances at a rate of around 90° 
per km. In addition, at ranges beyond about 2 km, the azimuthal component 0<j> B 

20 consistently lags the radial component 0 P B by around 25°. In the hydrocarbon-layer 
subterranean strata configuration, however, the behaviour is somewhat different. 
Beyond around 5 km, the azimuthal component 0«> R again lags the radial component 
0 P R by around 25°, however the phase of both components advances at a rate of around 
only 10° per km. This is significantly lower than that seen with the background 

25 subterranean strata configuration. Furthermore, at ranges between around 2 km and 5 
km, the difference in phase between azimuthal component 0 O B and the radial 
component 0 P R varies significantly. A phase separation anomaly is seen which varies 
from close to 0° phase difference between the radial and azimuthal components of 
detected electric field to a maximum of almost 60°. 

30 Figure 5B shows a graph plotting, for both the radial and azimuthal 

components, the difference in phase 0 R -0 B between the hydrocarbon-layer and 
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background subterranean strata configurations as a function of range R. The difference 
in the radial components is marked 9 P R -9 P B , and the difference in azimuthal 
components is marked 9 o R -0o R . The different rates of phase advancement seen with 
the hydrocarbon-layer and background subterranean strata configurations is apparent 
5 in the negative gradient of the curves beyond around 3 km. The differing behaviour at 
mid ranges (between around 2 km and 5 km) is apparent from the separation of the 
curves over this range. 

These differences in phase behaviour, namely the relatively slow advancement 
in phase of both radial and azimuthal components when a reservoir is present, and the 

10 strong range-limited variation in phase between the radial and azimuthal components 
seen at mid-ranges when the reservoir is present, provide two useful characteristics 
with which to determine the presence or absence of a hydrocarbon layer within an 
otherwise uniform background. 

For the practical application of controlled source electromagnetic surveying to 

15 hydrocarbon exploration, it is necessary that other common subterranean strata 
configurations do not lead to a behaviour similar to that seen in the hydrocarbon-layer 
subterranean strata configuration model. In particular, it is important to be able to 
distinguish between subterranean strata configurations which include a thin 
hydrocarbon layer and non-hydrocarbon containing subterranean strata configurations 

20 that have increasing resistivity with depth. 

Figure 6 shows in vertical section a highly schematic model of a non- 
hydrocarbon containing subterranean strata configuration. This subterranean strata 
configuration exhibits increasing resistivity with depth, which is a common feature of 
many large scale sedimentary structures. Due, for example, to increasing expulsion of 

25 conducting seawater with depth from rising overburden pressure. As with the 
background and hydrocarbon-layer subterranean strata configurations described 
above, in the non-hydrocarbon bearing subterranean strata configuration a section of 
seafloor 106 lies beneath a 10 km depth of seawater 104. The strata beneath the 
seafloor 106 comprise a series of sedimentary layers of increasing resistivity. A first 

30 layer 110 has a uniform resistivity of 1 Sim and a thickness of 1 km. A second layer 
1 13 has a uniform resistivity of 5 Qm and a thickness of 1 km. A third layer 1 14 has a 
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uniform resistivity of 50 £2m and a thickness of 1 km. Beneath the third layer 114 is a 
fourth layer 116 which has a resistivity of 100 Qm and extends downwardly for an 

infinite extent. A HED antenna 21 and a detector 25 are again shown as in Figure 3. 

l 

Figure 7 shows a graph which is similar to and will be understood from the 
5 description of Figure 5 A above, but which also includes modelled curves determined 
for the non-hydrocarbon bearing subterranean strata configuration. The modelled 
curves marked G P B , 8o B , 9 P R and 0<& R are the same as those shown in Figure 5A. The 
curves marked 9 P S and 9<> s show the radial and azimuthal components of the detected 
electromagnetic field seen with the non-hydrocarbon bearing subterranean strata 

10 configuration. 

The behaviour of the variation in phase of the detected electric field as a 
function of range beyond around 5 km is broadly similar for the hydrocarbon-layer 
subterranean strata configuration and the non-hydrocarbon bearing subterranean strata 
configuration. There are moderate differences in gradient between the models for the 

15 example shown, and in some circumstances this may allow the two subterranean strata 
configurations to be distinguished. (In other examples, the gradients are almost the 
same.) However, even if measurable, the value of the gradient is likely to be a fairly 
unreliable indicator of subterranean strata configuration in practice. This is because 
different absolute values of resistivity, for instance a more or less resistive 

20 . hydrocarbon layer in the hydrocarbon-layer subterranean strata configuration, or a 
more rapidly increasing resistivity with depth in a non-hydrocarbon bearing 
subterranean strata configuration, are likely to lead to changes in the observed 
gradients and cause confusion between the two models. 

However, at ranges between around 2 km and 5 km, there is nothing in the 

25 phase of the detected electric fields in response to the non-hydrocarbon bearing 
subterranean strata configuration which resembles the phase separation anomaly seen 
with the hydrocarbon-layer subterranean strata configuration. The phase behaviour 
seen with the non-hydrocarbon bearing subterranean strata configuration much more 
closely resembles that of the background subterranean strata configuration across this 

30 range. Accordingly, it is the phase separation anomaly (and not the gradient) which 
provides the most appropriate indicator of subterranean strata configuration. 
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Figure 8 shows a graph plotting the difference in the phase AG between the 
radial and azimuthal components of the detected electric field for the three model 
subterranean strata configurations described above as function of range R. The curve 
marked A9 B in Figure 8 represents the difference between the curves marked 0d> B and 
5 0 P B in Figure 7 (with negative values corresponding to the azimuthal component 
lagging the radial component). The curves marked A0 R and A0 S in Figure 8 
correspondingly represent the differences between the curves marked 0<d R and 0 P R , and 
0<d S and 9 P S in Figure 7 respectively. 

The phase separation anomaly seen with the hydrocarbon-layer subterranean 

10 strata configuration (curve marked A9 R in Figure 8) is apparent as a clear trough 
centred at a range of around 3.5 km. The magnitude of the phase separation anomaly 
at this point is almost 60°. This difference in phase between the radial and azimuthal 
components of the detected electric field is about 30° more negative than the largest 
difference seen with either the background or the non-hydrocarbon bearing 

15 subterranean strata configurations. With current technology, a phase differences 
between the radial and azimuthal components of the detected electric field of 10° can 
be clearly resolved. Accordingly, the presence or not of a trough similar to that seen in 
Figure 8 is easily detectable, and able to distinguish between a hydrocarbon-layer 
subterranean strata configuration of the type shown in Figure 3, and the model 

20 subterranean strata configurations shown in Figures 4 and 6. 

By distributing a linear array of detectors along a section of seafloor, and at 
each one recording suitable raw data in response to a source electromagnetic signal 
broadcast by a horizontal electromagnetic dipole source, plots such as those shown in 
Figure 8 can be generated from the phase information obtained from the raw data. The 

25 results of these plots can then be used to indicate the type subterranean strata 
configuration beneath a line joining the source and the detectors. Unlike previous 
survey methods, this can be done using a single dipole source and without the need 
to collect multiple data sets corresponding to different orientations of the source. 

By distributing a planar array of detectors on a section of seafloor, and at each 

30 one recording raw data in response to a source electromagnetic signal broadcast by a 
horizontal electromagnetic dipole source, plots such as those shown in Figure 8 can be 
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generated for a number of different directions once the phase information has been 
extracted from the raw data. Because of the insensitivity of phase to the azimuthal 
position of a detector with respect to the source dipole axis, the plots along each of the 
different directions achievable with a planar array of detectors can be obtained 
5 simultaneously, irrespective of the dipole source orientation. This allows a thorough 
two- or three-dimensional survey to be performed without even having to move the 
source. This contrasts to previous methods where a relatively long and complicated 
tow of the dipole source is required to utilise all of the detectors in a planar array, and 
then only with relatively low spatial sampling. Whilst all of the detectors in a planar 

10 array can be utilised without moving the source, in practical surveys employing the 
new method, it is likely that the source will nonetheless be moved, such as shown in 
Figure 2A. Each new source position provides an entire set of useful source-detector 
geometries, and so provides more comprehensive sampling of the subterranean strata 
configuration for a given number of detectors. In addition, by moving the source, 

15 surveys can be fully performed where the detectors are deployed over an area with a 
characteristic scale larger than the range of distances over which phase measurements 
can be reliably used to indicate the presence of a hydrocarbon layer. 

Figure 9A shows in plan view an area of seafloor 6 to be surveyed and which 
is similar to that shown in the prior art Figure L There are sixteen detectors 25 for 

20 recording the phase components described above. The detector are laid out in a square 
grid above a subterranean reservoir 56. Other detector distributions could be used 
instead, such as other grid shapes, or distributions that are not in a simple grid. (The 
constraints on detector placement patterns imposed by the amplitude-based methods 
of both Statoil and Sinha are therefore lifted) The subterranean reservoir 56 has a 

25 boundary indicated by a heavy line 58. The orientation of the subterranean reservoir is 
indicated by the cardinal compass points (marked N, E, S and W for North, East, 
South and West respectively) indicated in the upper right of the figure. To perform a 
survey using an embodiment of the new method, a source starts from location 'A' and 
is towed along a path indicated by the broken line 120 to location *B\ which marks 

30 the end of the survey path. At most points along the tow path, useful data can be 
collected from all of the detectors. For example, when the source antenna is at the 
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location marked by the reference numeral 57 in Figure 9A, all sixteen of the detectors 
25 are able to collect reliable data. This contrasts to the correspondingly similar 
location shown in Figure 1, again marked by the reference numeral 57, at which point 
no useful data can be collected using previous methods. Accordingly, the tow path 
5 shown in Figure 9 A, which is startlingly simple compared to that shown in Figure 1, 
actually provides a much greater amount of valid data. As noted above, with the new 
method, it is only necessary to know the relative positions of the source and detectors, 
and the orientation of each detector antenna such that the radial and azimuthal 
components of the detected electric field can be geometrically resolved. Since the 

10 orientation of the antenna is not critical, there is no need for the tow path 120 shown 
in Figure 9A to closely follow a pattern defined by the grid of detectors. In fact, it is 
preferable for the tow path to not align too closely with the north-south and east-west 
based detector grid, since for detectors in an end-on position (i.e. at an azimuth of 0° 
in the coordinate system shown in Figure 2B), the amplitude of the azimuthal 

15 component of the detected electric field will be small for a dipole source, and the 
phase of this component more difficult to accurately establish. In source-detector 
orientations where either the radial or azimuthal components of the detected electric 
field are small, other components of the detected electric field may be employed, for 
instance as described further below. 

20 Figure 9B is a graphical representation comparing the signal spread of the new 

phase method to the old inline/broadside amplitude method. The example reservoir 56 
bounded by the perimeter 58 is shown. The dipole source is at an arbitrary location 57 
within the reservoir with the HED antenna axis aligned W-E. In the old method, good 
quality inline amplitude data is only collectable within a narrow angular range 64 

25 indicated by W-E dark shading in the figure, and good quality broadside amplitude 
data is only collectable within a narrow angular range 62 indicated by the N-S dark 
shading in the figure. The angular ranges 62 and 64 need to be narrow to ensure that 
one of the inductive and galvanic signal components dominates over the other. Data 
collected by detectors in the four main quadrants 66 is essentially bad data to be 

30 rejected from analysis. On the other hand, in the new method, the situation is reversed. 
The broad quadrants 66 become the regions over which good quality data is collected, 



WO 2004/049008 



PCT/GB2003/005094 



-22- 

since they are the regions in which phase can be fully decomposed into the radial and 
azimuthal signal components needed for the phase-difference anomaly measurement, 
whereas the dark areas 62, 64 are angular areas where the collected data becomes 
unreliable since the magnitude of one of the radial and azimuthal signal components is 

5 likely to become too small causing signal-to-noise problems. 

The modelled phase responses shown in Figures 5A, 5B, 7 and 8 were all 
calculated for a horizontal electromagnetic dipole source transmitting a source 
electromagnetic signal at a frequency of 0.5 Hz. 

Figure 10 shows a graph plotting the modelled phase 9 of the radial and 

10 azimuthal components of the detected electric field for the hydrocarbon-layer 
subterranean strata configuration model shown in Figures 3 as a function of R for two 
different frequencies of source electromagnetic signal. The modelled radial and 
azimuthal components of the detected electric field seen in response to a dipole source 
transmitting at a frequency of 2 Hz are marked 0 p 2Hz and 9<j> 2Hz respectively, and the 

15 modelled radial and azimuthal components of the detected electric field seen in 
response to a dipole source transmitting at a frequency of 0.25 Hz are marked 9 p a25Hz 
and 6<j> 0,25Hz respectively. These curves, and also comparison with the modelled radial 
and azimuthal components of the detected electric field seen in response to a dipole 
source transmitting at a frequency of 0.5 Hz, marked 0 P R and 9<j> R in Figure 5A, 

20 indicate a frequency dependence to the characteristics of the phase separation anomaly 
indicative of a buried hydrocarbon layer. Towards higher frequencies, the phases of 
the radial and azimuthal components advance faster than at lower frequencies, and the 
scale over which the phase separation anomaly characteristic of a hydrocarbon layer's 
presence occurs is also seen to be frequency dependent. 

25 Figure 11A shows a graph plotting the difference in the phase A9 between the 

radial and azimuthal components of the detected electric field for the reservoir and 
background model subterranean strata configurations described above as function of 
range R, in response to a dipole source transmitting at a frequency of 2 Hz. The curve 

2Hz 

marked A9 R in Figure 11A represents the difference between the curves marked 9<p 
30 and 9 p 2Hz in Figure 10 (with positive values corresponding to the radial component 
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lagging the azimuthal component). The curve marked A0 B represents the 
corresponding data for the background model subterranean strata configuration. 

Figure 11B is shows a graph plotting the difference in the phase AG between 
the radial and azimuthal components of the detected electric field for the reservoir and 

5 background model subterranean strata configurations described above as function of 
range R, in response to a dipole source transmitting at a frequency of 0.25 Hz. The 
curve marked AG R in Figure 11B represents the difference between the curves marked 
g^o.25Hz and 0p O.25Hz - fl pi gure (with positive values corresponding to the radial 
component lagging the azimuthal component). The curve marked A0 B represents the 

10 corresponding data for the background model subterranean strata configuration. 

It is clear from Figures 8, 11A and 11B, that the range over which the phase 
separation anomaly occurs is smaller at higher frequencies. At 2 Hz (see Figure 11 A), 
the phase separation anomaly is centred at a range of around 3 km and occurs over a 
characteristic range of about 1 km. At 0.5 Hz (see Figure 8), the phase separation 

15 anomaly is centred at a range of around 3 km and occurs over a characteristic range of 
about 2 km. At 0.25 Hz (see Figure 11B), the phase separation anomaly is centred at a 
range of around 3 km and occurs over a characteristic range of about 3 km. For all 
frequencies, the maximum phase separation seen with the hydrocarbon-layer 
subterranean strata configuration is about 30° more negative than the phase difference 

20 that would be seen if the hydrocarbon layer were not present. This indicates that the 
presence of a hydrocarbon layer can be detected using a range of frequencies, each of 
which acts a probe of the subterranean strata configuration operating over a slightly 
different spatial scale. 

Figure 12A is a graph showing the effect of differing overburden thicknesses. 

25 The graph plots the difference in the phase AO between the radial and azimuthal 
components of the detected electric field for several hydrocarbon-layer subterranean 
strata configurations with different overburden thicknesses as function of range R. In 
this example, the source electromagnetic signal is at a frequency 0.5 Hz. Curves are 
plotted for different hydrocarbon-layer subterranean strata configurations which, while 

30 otherwise similar to the hydrocarbon-layer subterranean strata configuration shown in 
Figure 4, have overburden thicknesses of 0.25 km, 0.5 km, 1.0 km, 1.5km and 2.5 km. 
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The curves corresponding to each different overburden thickness are correspondingly 
marked in the figure. The curve marked 1.0 in Figure 12A is identical to the curve 
marked A0 R in Figure 8 since the overburden thickness in the model shown in Figure 
3 (and used for the modelling shown in Figure 8) is 1.0 km. The curve marked AG 3 in 
5 Figure 12A is similar to and will be understood from the similarly marked curve in 
Figure 8. For each of the curves corresponding to different overburden thicknesses, 
the magnitude of the phase separation anomaly is roughly similar, varying from about 
55° with an overburden thickness of 0.25 km to about 65° with an overburden 
thickness of 2.5 km. Accordingly, the method is equally able to detect a thin 

10 hydrocarbon layer at a range of depths beneath the sea floor. It is also apparent that the 
range at which the phase separation anomaly is maximum increases with increasing 
overburden thickness. This sensitivity of the range of maximum phase separation 
anomaly to overburden thickness can allow the depth of a reservoir to be determined 
with appropriate inversion modelling and suitable data coverage. 

15 Figures 12B, 12C and 12D are similar to and will be understood from Figure 

12 A. However, Figures 12B, 12C and 12D show the response of different overburden 
thicknesses to different frequencies of source electromagnetic signal. Figure 12B 
shows the response to a source electromagnetic signal at a frequency 0.25 Hz, Figure 
12C shows the response to a source electromagnetic signal at a frequency 1.0 Hz, 

20 Figure 12B shows the response to a source electromagnetic signal at a frequency 2 Hz. 
It can be seen that the phase separation anomaly is detectable with a range of 
frequencies over a range of overburden thicknesses. The magnitude of the phase 
separation anomaly is broadly similar at each of the different frequencies shown. As 
seen previously, the range over which the phase separation is apparent narrows with 

25 increasing frequency. 

Figure 13 A is a graph showing the effect of differing background resistivity. 
The graph plots the modelled phase 6 of the radial and azimuthal components of the 
detected electric field for a background subterranean strata configuration similar to 
that shown in Figure 3, but with different resistivity values for the uniform 

30 subterranean strata, as a function of range R. In this example, the source 
electromagnetic signal is at a frequency 0.5 Hz. The phase of the radial and azimuthal 
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components of the detected electric field are calculated for resistivity values of 1 Qm, 
5 Qm, 15 Qm, 50 Qm and 200 Qm, as marked in the figure. For each resistivity value 
the phase of the radial component of the detected electromagnetic field is shown as a 
dashed line, and the phase of the azimuthal component is shown as a solid line. The 
5 pair of curves corresponding to the 1 Qm resistivity value are identical to the curves 
marked 9 P B and 0o B in Figure 5A. It is clear that the resistivity value for the 
subterranean strata in a uniform model containing no hydrocarbon reservoir has a 
significant effect on the detected phase. The rate of advancement of phase with range, 
for both the radial and azimuthal components of the detected electric field, falls with 

10 increasing background resistivity. For example, at a resistivity of around 15 Qm, the 
rate of phase advancement with range is about 15° per km. This is similar to the rate 
of phase advancement seen with the reservoir subterranean strata configuration model 
and plotted in Figure 5A for ranges beyond around 5 km. This again demonstrates 
how absolute values of phase for each of the radial and azimuthal components can be 

15 an unreliable indicator of the likely presence of a hydrocarbon layer within an 
otherwise uniform resistivity background. 

Figure 13B shows a graph plotting the difference in phase AG between the 
radial and azimuthal components of the detected electric field for each of the different 
resistivity value background subterranean strata configurations shown in Figure 13A. 

20 Each curve is appropriately marked according to the resistivity value of the model to 
which it corresponds. The curve marked 1 Qm is identical to the curve marked A0 B in 
Figure 8. While the characteristic difference in phase for the radial and azimuthal 
components is depends on the resistivity of the uniform subterranean strata, none of 
the curves shown in Figure 13B display a range limited phase separation anomaly 

25 which, as seen in Figure 8, is indicative of the presence of a buried hydrocarbon layer. 
It is the resistivity contrast between a buried hydrocarbon layer and an otherwise 
uniform background which gives rise to the phase separation anomaly. Accordingly, 
by forming the phase difference between the radial and azimuthal components of the 
detected electric field in the manner described above, a hydrocarbon-layer containing 

30 subterranean strata configuration remains clearly distinguishable from a range of 
uniform subterranean strata configurations of differing resistivities. 
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Figure 14 is a graph showing the effect of hydrocarbon-layer resistivities. The 
graph plots the largest difference in phase max(A0) between the radial and azimuthal 
components of the detected electric field as a function of differing hydrocarbon-layer 
resistivity P, in a hydrocarbon-layer subterranean strata configuration which is 
5 otherwise similar to that shown in Figure 3. Curves are shown for source 
electromagnetic signal frequencies of 0.25 Hz, 0.5 Hz, 1 Hz and 2 Hz, as indicated in 
the figure. For example, with a source electromagnetic signal frequency of 0.5 Hz, and 
a hydrocarbon-layer resistivity of 100 flm, the largest difference in phase between the 
radial and azimuthal components of the detected electric field is about -58°. This 

10 particular value corresponds to the minimum seen in Figure 8 for the curve marked 
A0 R . Typical hydrocarbon-layer resistivities are between a few tens of Qm and a few 
hundreds of Qm. It can be seen from Figure 14 that for all typical values of 
hydrocarbon-layer resistivity, a difference in phase between the radial and azimuthal 
components of the detected electric field of at least 30° is seen for all source 

15 electromagnetic signal frequencies. At lower frequencies, it is even higher. 

This demonstrates that the above described method is able to detect 
hydrocarbon layers with different resistivities, and using a range of source 
electromagnetic signal frequencies. 

It can also be seen from Figure 14 that the maximum phase difference between 

20 the radial and azimuthal components of the detected field is greatest for a hydrocarbon 
layer with a resistivity of about 50 Qm. At resistivities above and below this value, a 
decreasing maximum phase difference is seen. 

It has thus been demonstrated that the presence of a hydrocarbon layer in a 
subterranean strata configuration can be detected by observing the phase difference 

25 between radial and azimuthal components of detected electric field in response to a 
source electromagnetic signal from a horizontal electric dipole source. This has been 
shown to work over a wide range of source frequencies, for differing depths of burial 
of a hydrocarbon layer and for different subterranean strata configuration resistivity 
values. 

30 

Alternative Embodiments 
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Whilst in the above examples the radial and azimuthal components have been 
considered, similar techniques can also be employed using different components of 
the detected electromagnetic field. For instance, if the detectors 25 shown in Figure 2 
5 were configured to also record the phase of the vertical component of the detected 
electric field (i.e. perpendicular to both the radial and azimuthal directions), the 
vertical component could be used in combination with another component to probe a 
subterranean strata configuration. 

Figure 15 A shows a graph plotting the modelled phase 0 of the vertical and 
10 azimuthal components of the detected electric field for both the background 
subterranean strata configuration and the hydrocarbon layer subterranean strata 
configuration models shown in Figures 3 and 4 respectively as a function of R. The 
phase is measured relative to a source electromagnetic signal transmitted by the HED 
antenna. In this example, the source electromagnetic signal is at a frequency 0.5 Hz. 
15 The vertical and azimuthal components of the detected electric field for the 
background subterranean strata configuration are marked 8 Z B and 0a> B respectively and 
the corresponding components of the detected electric field for the hydrocarbon-layer 
subterranean strata configuration are marked G Z R and 0<d R respectively. Except for 
showing the vertical rather than the radial components of the detected electric field, 
20 this figure directly corresponds to Figure 5A. 

Figure 15B shows a graph plotting the difference in the phase A0 between the 
vertical and azimuthal components of the detected electric field for the two model 
subterranean strata configurations included in Figure 15 A. The curve marked A0 B in 
Figure 15B represents the difference between the curves marked 0<p B and 0 Z B in Figure 
25 15 A (with negative values corresponding to the vertical component lagging the 
azimuthal component). The curve marked A0 R Figure 15B correspondingly represents 
the difference between the curves marked 0<j> R and 0 Z R in Figure 15 A. 

Figure 16A shows a graph plotting the modelled phase 0 of the vertical and 
radial components of the detected electric field for both the background subterranean 
30 strata configuration and the hydrocarbon-layer subterranean strata configuration 
models shown in Figures 3 and 4 respectively. The phase is measured relative to a 
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source electromagnetic signal transmitted by the HED antenna. In this example, the 
source electromagnetic signal is at a frequency 0.5 Hz. The vertical and radial 
components of the detected electric field for the background subterranean strata 
configuration are marked 0 Z B and 8 P B respectively and the corresponding components 
5 of the detected electric field for the hydrocarbon-layer subterranean strata 
configuration are marked 0 Z R and 6 P R respectively. Except for showing the vertical 
rather than the azimuthal components of the detected electric field, this figure directly 
corresponds to Figure 5 A. 

Figure 16B shows a graph plotting the difference in the phase A9 between the 

10 vertical and radial components of the detected electric field for the two model 
subterranean strata configurations included in Figure 16A. The curve marked A0 B in 
Figure 16B represents the difference between the curves marked 0 P B and 9 Z B in Figure 
16A (with negative values corresponding to the vertical component lagging the radial 
component). The curve marked A8 R Figure 18B correspondingly represents the 

15 difference between the curves marked 9 P R and 0 Z R in Figure 18 A. 

Figures 15B and 16B both indicate that the difference in phase between the 
vertical component of detected electric field and either of the azimuthal or radial 
components is also sensitive to the presence of a hydrocarbon layer in an otherwise 
uniform background subterranean strata configuration. The phase difference seen 

20 between the vertical and azimuthal components displays both a negative and a 
positive lobe compared to the background subterranean strata configuration with a 
crossover at a range of around 4 km. This would be particularly useful indicator for 
use in survey areas where background resistivity is poorly constrained. The qualitative 
behaviour of the phase of the vertical component of the detected electric field is 

25 approximately similar to that of the azimuthal component. However, a larger phase 
separation is seen when comparing the radial component with the azimuthal 
component than when comparing the radial component with the vertical component. 
Accordingly, the azimuthal component will generally be preferred when comparison is 
made with the radial component, unless, for instance, the magnitude of the azimuthal 

30 component is small, for example, where a detector is very close to an end-on 
orientation. 
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In the above description, and in Figures 5A, 7, 10, 13 A, 15A and 15B, the 
absolute phase of various components of the detected electric field has been 
considered relative to the source electromagnetic signal phase. However, in practice, 
since it is the relative between different components of the detected electric field seen 
5 at the detector which is indicative of the presence of a hydrocarbon layer, the detected 
components may be directly compared without reference to the phase of the source 
electromagnetic signal. 

Finally it will be understood that the invention is equally applicable to 
surveying of freshwater, for example large lakes, so that references to seafloor, 
10 seawater etc. should not be regarded as limiting. 

Summary 

15 It has been demonstrated how a phase separation anomaly occurs in response 

to a hydrocarbon layer which is not seen with a background subterranean strata 
configuration. This allows the detection of subterranean hydrocarbon reservoirs and 
hydrocarbon bearing layers. The technique has many advantages over previous 
methods, for example: 

20 • Previous techniques based on comparison of amplitude measurements require 
the collection of both end-on and broadside data for each receiver to be reliable. This 
requires multiple orthogonal survey tow paths (see Figure 1). Using a technique such 
as described above, a survey may be completed more thoroughly with a much shorter 
and less complex tow path (see Figure 9 A). 

25 • The orthogonal towpaths required by previous methods lead to sampling of 
different parts of a target structure. Because only single source-receiver pairs are 
required for a phase-based detection of the reservoir, there is reduced interpretational 
ambiguity arising from the dimensionality of the target structure. This overcomes the 
limitation of the prior art Sinha method in which in-line data from a given receiver 

30 will come from one source position and the corresponding broadside data at the same 
range will generally come from a different source location. This means that the 
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stmcture sampled between the source and detector for the compared in-line and 
broadside data will not be the same. With the new method, this problem does not 
arise, since all the data is collected from a single source position, so both phase 
components in the processed signal are derived from sampling the same structure. 
5 • The above described technique is almost independent of detector azimuth 
relative to a source's dipole axis. Since the method is less dependent on the 
orientation of the source, geometry-related errors are much reduced. In order to 
decompose the detector signals into radial and azimuthal (or whichever components 
are desired) it is only necessary to know the relative positions of the source and 
10 detectors, and the orientation of the detector antenna. These can be easily determined 
using existing technology. 

• The phase separation seen above is range limited and can be controlled by 
varying the frequency of the electromagnetic source. If the source were to broadcast at 
several discrete frequencies (either by employing multiple source antenna or a tuned 

15 source for example) improved vertical resolution can be achieved. 

• For a particular source dipole transmission frequency, the range dependence of 
phase separation can be used to indicate the depth to the resistive layer. 

• Phase data are relatively insensitive to structures which are local to the 
receiver. 
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